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ABSTRACT 

We discuss the importance of feedback via photoionization and Compton heating on 
the co-evolution of massive black holes (MBHs) at the center of spheroidal galaxies 
and their stellar and gaseous components. We first assess the energetics of the ra- 
diative feedback from a typical quasar on the ambient interstellar medium. We then 
demonstrate that the observed Mbh-c relation could be established following the con- 
version of most of the gas of an elliptical progenitor into stars, specifically when the 
gas-to-stars mass ratio in the central regions has dropped to a low level ~0.01 or less, 
so that gas cooling is no longer able to keep up with the radiative heating by the 
growing central MBH. A considerable amount of the remaining gas will be expelled 
and both MBH accretion and star formation will proceed at significantly reduced rates 
thereafter, in agreement with observations of present day ellipticals. We find further 
support for this scenario by evolving over an equivalent Hubble time a simple, physi- 
cally based toy model that additionally takes into account the mass and energy return 
for the spheroid evolving stellar population, a physical ingredient often neglected in 
similar approaches. 
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1 INTRODUCTION 

Elliptical galaxies are, with few exceptions, poor w ith re- 
spect to interstellar gas (e.g. O'Sulliv an et alj l200lf) . Also, 
elliptical galaxies invariably contain central massive black 
holes (MBHs) , and there exists a tight relationship between 
the characteristi c stellar velocity dispersion a and the black 
hole mass Mbh iFerrarese fc MerrittlfeOOct iTremaine et alJ 
120021) . and between Mbh and th e host spheroid mass in 
stars, M, llMagorrian et al.lll998l) . Are these two facts re- 
lated? Here we focus on a scenario in which the mass of the 
central MBH grows within gas rich elliptical progenitors un- 
til star formation has reduced the gas fraction in the central 
regions to of order 1 per cent of the stellar mass. Then ra- 
diative feedback, during episodes when the luminosity from 
the central MBH approaches its Eddington limit, heats and 
drives much of the remaining gas out of the galaxy, limiting 
both future growth of the MBH and future star formation 
to low levels. 

We show that for a typical quasar spectral energy dis- 
tribution (SED) the limit on the central MBH produced 



by this argument coincides accurately with the observed 
Mbh-t relation and we predict that the observed power 
law should break down below 10 4 Mq and above a few 
1O 9 M . In our calculatio ns, w e adop t the average quasar 
SED derived bv lSazonov et alJ J2004L hereafter SOS) from 
cosmic background radiation fields supplemented by in- 
formation from individual objects. Of key importance is 
that the UV and high energy radiation from a typical 
quasar can photoionize and heat a low density gas up to 
an equilibrium Compton temperature (Tc ~ 2 x 10 7 K) 
that exceeds the virial temperatures of giant ellipticals; the 
radiative effects on cluster gas are expected to be min- 
imal (e.g. ICiotti. Ostriker fc Pellegrini! 12004) but the ef- 
fects on gas within the host galaxy can be substantial. We 
note that the pre s ent paper i s complementary to those by 
ICiotti fc Ostrikerl <1997t 1200 it hereafter CO97,01) in that, 
while it does not attempt to model the complex flaring be- 
havior of an accreting MBH with an efficient hydrodynami- 
cal code, it does do a far more precise job of specifying the 
input spectrum and the detailed atomic physics required to 
understand the interaction between that spectrum and the 
ambient interstellar gas in elliptical galaxies. 
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The evolutionary scenario considered could explain sev- 
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oral key observational facts. First, that giant ellipticals are 
old - they end their period of star formation early in cos- 
mic time, since the radiative output from the central MBHs 
limits (in cooperation with the energetic input due to star 
formation) the gas content to be at levels for which ongo- 
ing star formation is minimal. Secondly, gas rich, actively 
star forming galaxies at redshift z ~ 3, including Lyman 
Break Galaxies and bright submillimeter SCUBA galax- 
ies, generally ex hibit only moderate active galactic nucleus 
(AGN) activity iSteidel et aljbooi lAlexander et alJl2003t 
iLehmer et all l2004t) . indicating that their central MBHs 
are still growing. This suggests that the formation of a 
spheroid probably closely preceeds a quasar shining phase, 
which finds further support in spectroscopic observations 
showi ng that quasars live in a metal-enriched environment 
(e.g. lHamann fc Ferlandl fl999h . The redshift evolution of 
the quasar emissivity and of the star formation history in 
spheroids is thus expected to be roughly paral lel since z ~ 3, 
which also is consistent with observations (e.g. lHaiman et al] 
l2004lHeckman et aljEool . 

We note that many authors have already recognized the 
importance of feedback as a key ingredie nt of the mutual 
MBH and gala xy evolution (amon g them [ B inney & Taborl 
19951: CQ97: ISilk fc Reel \l99& iFabianl Il999t CO01; 
' Burkert fc Silkl 1200 J ICiotti fc van Albadal 12001 



TT 



Gran ato et alJl2Q0ll:ICavaliere fc Vittorinl20 02: Kine 2003; 
Wvithe fc Loebl 120031: lHaiman et alJ 12004 iGranato et al ' 



2004 iMurrav et al.ll2004 ISpringel. Pi Matteo fc Hernauisl 
2004). What is new about this work is the stress on one 
component of the problem that has had relatively little 
attention - the radiative output of the central MBH. 

We emphasize that the radiative output is not the only, 
nor even necessarily the dominant mechanism whereby feed- 
back from accre tion onto central MBHs can heat gas in el- 
liptical galaxies. iBinnev fc Taborl 1119951) have stressed that 
the mechanical input from radio jets will also provide a sig- 
nificant source of energy, and much detailed work has been 
performed to follow up this suggestion. Strong evidence has 
accumulated over recent years that the radiative losses of 
hot gas in the 100 kpc-scale cores of galaxy clusters may be 
counteracted by the mechanical input from the central gi- 
ant ellipti cal galaxies, which usu ally host a low luminosity 
AGN (e.g. IChurazov et aljl20o3) . It is however not obvious 
that all AGNs, and in particular luminous quasars, produce 
radio jets, whereas all do appear to produce high energy 
radiation. In any case, the two mechanisms are complemen- 
tary, and in this paper we are exploring only the radiative 
feedback. Given what we now know about the prevalence 
of MBHs in ellipticals, the tight relationship between the 
MBH mass and stellar mass, and the known efficiency for 
converting accreted mass into electromagnetic output with 
a specific spectrum, the consequence of radiative feedback 
can now be estimated with reasonable accuracy, more easily 
than can be acheived for the mechanical energy feedback. 

This paper is organized as follows. In Section ^. ll we in- 
vestigate the possibility of radiatively preheating a spherical 
accretion flow within 1-100 pc from a MBH. In Section [2.21 
we assess the conditions for a significant heating of inter- 
stellar gas at larger distances from the central quasar. In 
SectionE]we propose that radiative feedback could be a key 
factor leading to the observed Mbh-t relation. In Section 2] 
we elaborate on this idea in the context of the MBH-galaxy 



co-evolution, using a simple, but physically motivated, one- 
zone model. 



2 ENERGETICS OF THE RADIATIVE 
FEEDBACK 

Radiative feedback from accretion onto a MBH can affect 
the accretion flow within the MBH sphere of influence, i.e. 
in the central 1-100 pc of the host galaxy [see equations (1141 
and 1161 1. as well as the interstellar gas at larger distances. 
Adopting the average quasar SEP from SOS, we consider the 
effects arising in the former and later case in Section ^. H and 
12.21 respectively. We note that significant feedback effects 
are expected only when the luminosity L of the central MBH 
approaches the Eddington luminosity 



L Edd = 1.3x10 Y^ergs , 

hence our use of a SEP typical of quasars throughout. 



(1) 



2.1 Preheating of a spherical accreting flow onto a 
MBH 

Using hydrodynamical simulations and fiducial quasar 
SEPs, CO01 demonstrated that accretion of gas onto the 
central MBHs in elliptical galaxies cannot proceed stably at 
sub- Eddington rates if the mass-to-radiation conversion effi- 
ciency e (= A£ 7 /Amc 2 ) is high (0.1), as is characteristic of 
quasars. Puring high luminosity episodes, the high energy 
radiation, primarily X-rays, from the MBH unbinds the ac- 
creting gas through Compton heating, which temporarily 
leads to a decrease in the mass inflow rate and consequently 
in the central luminosity. The SEPs utilized by CO01 cov- 
ered a fair range, in general corresponding to a high Comp- 
ton temperature. Our aim here is to verify whether the con- 
clusions remain unaltered for the average quasar SEP once 
photoionization is fully taken into account. We note that our 
treatment below is similar to previous semi-analytic stud- 
ies of preheating of quasi-spherical accreting gas flows (e.g. 
[ Cowie et al]ll97St iKrolik fc Londonlll983t iPark fc Ostrikerl 

SOS (see their fig. 5) computed for the average quasar 
SEP the equilibrium temperature T eq of gas of cosmic chem- 
ical composition as a function of the "ionization" parameter 



L 



(2) 



where n is the hydrogen number density and r is the distance 
from the MBH (note the use of the bolometric luminosity in 
the definition of £). In the temperature range 2 x 10 4 -10 7 K, 
to a good accuracy 



T cq (C) » 2 x 10 2 £K, 



(3) 



while for £ < 100 and £ > 5 x 10 4 , T oq « 10 4 K and 
2 x 10 7 K, respectively. T cq is the temperature at which 
heating through Compton scattering and photoionization 
balances Compton cooling and cooling due to continuum 
and line emission, on the assumption that the plasma is in 
ionization equilibrium, appropriate for the problem under 
consideration. In Fig.Qwe plot lines of constant T eq on the 
(r, nLvdd/L) plane for Mbh = 10 8 M© and Mbh = 1O 9 M 
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[note that the full expression of T e q(£) is used for this pur- 
pose]. 

Suppose now that the gravitational potential is due to 
the central MBH alone. Then the condition 



5 fcTcq(?) _ GAWmp > Q 



(4) 



(here ll = 0.61 is the mean molecular weight) defines a sit- 
uation where gas of density n, located at r, will be heated 
to T eq by the central radiation and blown out of the MBH 
potential . For given Mbh, L/L^dd and r, gas with density 
below a certain critical value, as implicitly given by equa- 
tion @, cannot accrete onto the MBH, which corresponds 
to the dashed areas in Fig. 

We now specify the general conditions above to Bondi 
accretion. Note that we are not assuming Bondi accretion 
onto the central MBH, since the infalling material is ex- 
pected to possess non-negligible angular momentum and 
form an accretion disk or a more complicated central flow 
near the MBH 2 . We thus regard the Bondi relations as an in- 
dication of the flow rates into the central regions as a source 
of mass for the central flows. 

The Bondi flow is transonic inside the radius 



Rb 



27/cT 



1fi 1 Af BH ( T \ 

PC 7 1O 8 M Vio 6 kJ 



(5) 



where T is the temperature at _Rb and the sound speed 
c s = yj ykT/ (im p . We shall assume 7 = 1 (isothermal flow) 
for estimations below, the results depending very weakly on 
7. Comparison of equation I^J with equation @ shows that 
Bondi accretion of gas at temperature T can be disrupted if 
the central luminosity L is sufficiently high that T c<1 (Rb) > 
T, or equivalently L > L cr it, where 



L clit (T) = i{T)Rl{T)n{R B ) 



(6) 



and £(T) is the ionization parameter corresponding to T cq = 
T. If L > L cr it, the gas in the vicinity of Rb will experience 
net radiative heating. On the other hand, the central lumi- 
nosity is determined by the pre-existing conditions at Rb 
(so far as L < Z/Edd): 



x 4TvRi(T)c s (T) x nm p n t (R E 



(7) 



where n t ~ 2.3n is the total particle number density. 

Fig. |21 (upper panel) shows the ratio L/L cr it as a func- 
tion of gas temperature for three values of the accretion 
efficiency, e = 0.01, 0.05 and 0.1. The two higher values 
would correspond to accr etion via a standard thin disk 
JShakura fc Sunvaevll973l) . The presented curves have been 
obtained from equations ||SJ and using the full expression 



1 We note that in the case of the MBH emitting at near L^dd, 
the radiation pressure will effectively reduce the MBH gravita- 
tional potential, making possible an outflow in the directions not 
screened by the central accretion disk for lower values of £ than 
implied by equation |IJ. 

2 In the case of quasars, it is believed that accretion occurs via 
a geometrically thin, optically thick disk, the presence of which 
naturally explains both the fact that most of the radiative output 
takes place at ultraviolet wavelengths and the inf erred high radia- 
tive efficiency (e f» 0.1) of accretion (e.g. lYu &: Tremaindl20o3) . 
We briefly discuss the consequences of the presence of a central 
accretion disk at the end of this Section. 




Figure 1. Energy-based gas escape condition near a black hole of 
mass Mbh in the (r, n) plane, under the assumption that the gas 
is in ionization equilibrium. Slanted lines correspond to constant 
values of the ionization parameter as given by equation (3) for dif- 
ferent equilibrium gas temperature T eq (labels along the lines). 
T eq is shown up to the Compton temperature (1.9 X 10 7 K) char- 
acterizing the average quasar. The dashed area below the solid 
curve given in equation J3J corresponds to gas configurations en- 
ergetically unbound. In the upper and lower panel Mbh = 10 s Mq 
and Mbh = 10 9 Mq, respectively. 



of T eq (£), and can be approximated in a limited temperature 
range by the expression 

-1/2 



L 



6- 



T 



2 x 10 4 K < T < 10 7 K, 



(8) 



Lcrit 0.1 U0 6 K, 

which follows from equation ©. Note that the ratio L/L cr it 
does not depend on Mbh- The fact that L > L cr it over most 
of the plot in Fig. [5] suggests that stationary accretion at 
high rates is impossible, unless the gas supplied into the 
Bondi sphere from outside is hotter than the Compton tem- 
perature of the average quasar SED, i.e. T > Tc ~ 2x 10 7 K. 
If T < Tc, the inflowing gas near R B will be heated and 
driven out by the quasar radiation, and the MBH luminosity 
will be reduced for the time required for significant amounts 
of gas to cool and accumulate again in the central regions, 
i.e. the result will be relaxation oscillations. 

However, such a dramatic heating effect as described 
above will only occur if the characteristic flow time at _Rb, 



^dyn 



Rb 



(2GMbk/RbY' 2 j 10 8 M© 

is longer than the heating time scale 
_ n t (R B )kT 

thcat = pj-y ^ • l-LUJ 

Here T(T,£) is the plasma net heating (cooling) rate per 
unit volume, which in the case of the average quasar SED 
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Figure 2. Top panel: the ratio L/L CI i t as a function of T for 
Bondi accretion, as given by equations JSJ and Q. The label 
above solid curves represents the adopted radiative efficiency. 
L/L crl t > 1 implies that stationary accretion is impossible due 
to strong preheating. Bottom panel: ratio of the characteristic 
heating and flow times at the Bondi radius, as given by equa- 
tions £3 and CSJ: 

as a function of the gas temperature T for 
e = 0.1 and different values of L/-LEdd (labels). Solid lines cor- 
respond to net heating, while dashed lines to Compton heating 
only. iheatAdyn > 1 implies that radiative heating is inefficient. 



can be estimated from formulae given in Appendix IA3.3I 
Assuming e = 0.1, we find that for L/ L cr n approximated by 
equation J§) 



r w 3 x 10~ 23 n 2 



(i?B)ergs 'cm 3 , 



the 



the 



2 x l(f yr 



T 



tdy 



0.03 



(l0 6 K 
L \ _1 T 



- 1 / 2 Mb 



(— y 

V Z/Edd / 



V I/Edd / 



10 B K 



1O 8 M 

2 x 10 4 K < T < 10 7 K. (11) 



A more accurate computation leads to the curves shown in 
Fig. |21 (lower panel), giving the ratio iheat/tdyn as a function 
of T for different values of L/L^dd- Also plotted is the ratio 
tc/tdyn for the case of pure Compton heating/cooling. One 
can see that Compton preheating is important when T > 
10 6 K and L £ 0.1L E dd- For T < 10 6 K, preheating due to 
photoionization is significant at even lower accretion rates. 

The above discussion suggests that if quasars are fed 
by warm (T < 10 7 K) interstellar gas from distances ~10- 
100 pc, their emission is expected to be significantly variable 
on time scales ~ td yn ~10 4 -10 5 years. Only for gas tem- 
peratures approaching 10 7 K, is stationary accretion at sub- 
Eddington rates possible (see Fig. 

We note however that, if the final stages of accretion 
occur via a geometrically thin or slim disk, as expected 
for quasars (see Footnote |5J , the above conclusion should 



change as regards the variability time scale. Indeed, the char- 
acteristic viscous time in a geometrically thin, opticall y thick 
accretion disk is given by jShakura fe Sunvaevlll973T) 



ivisc(r) ~3x 10yr 



\0.lj \L E dd) \10 s Mq J 



0.01 pc 



(12) 



where a is the viscosity parameter. The above formula is 
valid in the innermost, radiation pressure dominated zone 
of the disk, where the main source of opacity is Thomson 
scattering. Such disks, characterized by a radially constant 
accretion rate, are unlikely to extend beyond 10~ 3 -10~ 2 pc 
because they would become s elf-gravitat i ng ifaiomrel Il964l : 
iKolvkhalov fe Sunvaevl ll98Ct iGoodmanl 120031) . Since the 
quasar life time is expected to exceed 10 7 yr, the accretion 
disk must be constantly replenished with gas of small angu- 
lar momentum (presumably from a Bondi- type flow). 

In the case of quasars (Mbh J> 1O 8 M0, L ;> 0.11/Edd), 
self-gravity is expected to truncate the disk already in the 
innermost zone, at the distance 



/ a \ 2 / 9 / L \ 4/9 / 
^.0.014 pc (_) ( ) 



M 



\^1O 8 M 



7/9 



(13) 



and the disk is expe cted to be stable against the therma l- 
viscous instability l|Burderi. King fc Szuszkiewiczl [l998). 
Comparison of equations © and 11211 implies that the 
timescale on which accreting matter drifts from the outer 
boundary of the disk to the MBH, t v i sc (R sg ), can exceed 
the characteristic dynamical time of the external Bondi-type 
flow, tdyn(^B). In such a case, variations in the gas inflow 
rate at Rb as a result of preheating will affect the MBH ac- 
cretion rate with a time delay ~ tvisc, leading to variations 
in the quasar luminosity also on time scales ~ tvisc- 

2.2 Heating of ISM in spheroids 

The gravitational potential of the central MBH is overcome 
by the potential of the host spheroid at distances larger than 
the MBH radius of influence 



Rbh 



GMbh 



in Mbh ( — V 
PC 10 8 M Q UOOkms-V 



(14) 



'1O 8 M 

where a is the characteristic (virial) one-dimensional stel- 
lar velocity dispersion in the galaxy. Comparison of equa- 
tion 11411 with equation © implies that Rb < Rbh when 
T £ T vir , where 



k 



3.0 x 10 B K 



200 km s" 



(15) 



is the galaxy virial temperature. Therefore, our discussion of 
preheating in i|2.1| pertains to situations where the accreting 
flow has a temperature comparable to or higher than the 
galaxy virial temperature. 

Below we assess the conditions required for the cen- 
tral MBH to significantly heat the interstellar gas over a 
substantial volume of the galaxy, regardless of the type of 
accretion flow established in the central regions. In this Sec- 
tion we shall assume that the MBH has a mass as given by 
the observed Mbh-c r elation for local ellipticals and bulges 
jTremaine ct al. 20q3): 
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M B h = 1.5 x l(fM, 



( °- V 

1 200 km s" 1 / 



(16) 



.200 km s- 1 . 

Note that this assumption will be dropped in SJ3J where we 
attempt to predict the Mbh-ct" relation. 

From equations (1151 and 11611 we can find the critical 
density n cr it, defined by 



T C q(L/n cr i t r ) — T v i r 



(17) 



as a function of distance r from the galaxy center. Gas with 
n < n CI it will be heated above T v i r and expelled from the 
galaxy, whereas gas cooling will dominate over Compton 
and photoionization heating if n > n C rit. We show in Fig. |3] 
the resulting heating diagrams on the (r, n) plane for two 
galaxies, with a — 180 and 320 km s -1 , corresponding to 
Mbh = 10 s and 10 9 M©, respectively. Also shown in Fig. 
is the n(r) line at which the cooling time tc»oi,o of gas at 
r v i r in the absence of radiative heating equals the charac- 
teristic dynamical time, td yn = r/a. Gas with density above 
this line will cool down before reaching the galactic nucleus 
if there are no heating mechanisms other than gravitational 



compression. Comparison of the T e , 



Tvir line with the 



tcooi.o — tdyn one implies that if the gas within a few kpc of 
the MBH is sufficiently tenuous to sustain a subsonic cool- 



ing flow (so that t c 



> idyn and T « Tvir), the central 



quasar emitting L > 0.1Z/Edd will be able to heat the gas 
above T v j r . 

In practice, provided that T cq > T v i r , significant heat- 
ing will take place only out to a certain distance that de- 
pends on the luminosity and duration of the quasar out- 
burst. Since the MBH releases via accretion a finite total 
amount of energy, eMbhc 2 , there is a characteristic limiting 
distance given by: 



Rc = ( 



(TtcMbh \ 1/2 



37rm c 

500pc (oT) 1/2 ( 



400pc (cfr) 



200 km s" 1 



1/2 



Mb 



10 s M, 



1/2 



(18) 



Inside this radius, each electron-proton pair will have re- 
ceived at least 3fcTc ~ 6keV of energy through Compton 
scattering of hard X-rays from the MBH when it has ac- 
creted mass Mbh- is defined assuming that the only 
heating mechanism is Compton scattering, and thus per- 
tains to the limit of low density gas, fully photoionized by 
the radiation from the MBH. 

More relevant for the problem at hand is the distance 
out to which low density gas will be Compton heated to 
T > T ■ ■ 



Ri = Rc 



1/2 



(19) 



200 km s" 1 ' 

Yet another characteristic radius is the one within which 
gas of critical density n cr i t will be heated to T > T V [ T by 
photoinization and Compton scattering: 

1/2 

(20) 



R 2 = Ri 



r(n c 



where Fc and T are the Compton and total heating rates, 
respectively. Depending on the gas density (0 < n < n cr i t ), 
the outer boundary of the "blowout region" will be located 
somewhere between Ri and i?2. The size of the heating zone 
may be compared (see Fig.[3Jl with the galaxy effective (half- 




r (pc) 



Figure 3. The (r,n) plane for a galaxy with a = 180 km s 1 
(T vir = 2.4 X 10 6 K, Af B H = 10 8 M Q , upper panel), and with 
a = 320 kms" 1 (T vir = 7.7x 10 6 K, M B h = 10 9 A/ e , lower panel). 
In the dashed area gas can be heated above T v i r by radiation 
from the central MBH emitting at the Eddington luminosity. The 
upper boundary of this area scales linearly with luminosity (as 
shown by the dashed line corresponding to L = O.lLEdd)- From 
left to right vertical lines correspond to i?BH [equation 1141 1 , Rc 
[equation 1181 1 . Ri [equation 1191 1 . R2 [equation 1201 1. and Rc 
[equation 1211 1. The thin slanted line bounds from the above the 
zone where the intrinsic cooling time of gas at T v ; r is longer than 
the dynamical time. See text for further explanation. 



light) radiu s R e . We estimate from the results of lFaber et alJ 
Jl997l) and lBernardi et"all i2003l) that for early- type galax- 



Re 



• 4,000 pc ( 



200 km s- 



(21) 



and note that in reality the correlation between R e and a is 
fairly loose. 

The different characteristic distances defined above are 
shown as a function of Mbh in Fig. 2] One can see that the 
radiative output of a black hole of mass < 10 7 Mq can un- 
bind the interstellar gas out to several R e in relatively low 
luminosity elliptical galaxies. In the case of more massive 
black holes/galaxies with Mbh ~1O 8 -1O 9 M0, the heating 
will be localized to the innermost 0.3-0.5ii c ; the gas fur- 
ther out can be heated only indirectly, presumably through 
shock waves propagating from the radiatively heated cen- 
tral regions (see CO01). Remarkably, R2 ~ 2-4 kpc, almost 
independently of A/bh- This results from a combination of 
two opposite trends: the lighter the MBH is, the less to- 
tal energy it emits, and the lower T V [ T is, the more im- 
portant becomes the role of photoionization heating com- 
pared to Compton heating for gas with n<n cr it- We point 
out that radiation characterized by the average quasar SED 
could not unbind the gas in galaxies with a > 500kms~ 
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Figure 4. Different characteristic heating radii, Rc (dotted line), 
Rl (short-dashed line), and R2 (long-dashed line), as a function 
of A/bH) compared with the general trend for the effective radius 
of elliptical galaxies (solid line). 

(Mbh > 6 x 1O 9 M0) if such giant galaxies existed, because 
then T vir > T c » 2 x 10 7 K. On the contrary, the gas in the 
central regions of such galaxies would be Compton cooled by 
the radiation from the central quasar. Presumably it is a co- 
incidence that the upper limit to the velocity dispersion for 
real elliptical galaxies is close to the Compton temperature 
found for the radiation from typical massive quasars. 

We finally note that in the case of a quasar outburst 
with given L and duration t, during which a relatively small 
amount of mass AMbh <C Mbh is accreted, the radiative 
heating front will propagate out to 



Ri, 2 (L,t) = Ri, 2 



1/2 



t 



2 x 10 7 yr 



1/2 



ai) <*> 



in time t, where R.1,2 are given by equations 1)19^ . (1201 . It is 
worth noting that in the case of Eddington- limited accretion, 
the characteristic heating time at R\,2 is 2 x 10 7 (0.1/e) yr, 
or approximately half the Salpeter time scale (time required 
for the black hole mass to double). 



3 THE PROPOSED ORIGIN OF THE M B n-a 
RELATION 

We now address the central issue of this work, namely the 
possibility that radiative feedback played the key role in 
establishing the observed Mbh-o" relation. 

Below we elaborate on the following general idea. Before 
the MBH grows to a certain critical mass, MBH,crit, its radi- 
ation will be unable to efficiently heat the ambient gas, and 
accretion onto the MBH will proceed at a high rate. Once the 
MBH has grown to Mbh, crit, its radiation will heat and expel 



a substantial amount of gas from the central regions of the 
galaxy 3 . Feeding of the MBH will then become self-regulated 
on the cooling time scale of the low density gas. Subsequent 
quasar activity will be characterized by a very small duty 
cycle (~0.001), as predicted by hydrodynami cal simulations 
(CQ97, 01) and suggested by observations llHaiman et alJ 
l2004lHeckman et alfeooll) . MBH growth will be essentially 
terminated. 

Suppose that the galaxy density distribution is that of 
a singular isothermal sphere, with the gas density following 
the total density profile: 



(r) = 



M e 



(23) 



At 2-nGr 2 ' 

Here M gas and At are the gas mass and and total mass 
within the region affected by radiative heating. The size of 
the latter is uncertain but is less than a few kpc (see H2.2f . 
so that At is expected to be dominated by stars (M*<M) 
rather than by dark matter. 

Radiation from the central MBH can heat the ambient 
gas up to the temperature 



..5xl0 3 Krr^- 



(M gaa y 
V At J 



Mb 



/ 200 km s" 



10 s M Q V 



■(24) 



This approximate relation is valid in the range 2 x 10 4 -10 7 K 
and follows from equations (|3J and 1231 . Remarkably, T eq 
does not depend on distance for the adopted r~ 2 density 
distribution. We then associate the transition from rapid 
MBH growth to slow, feedback limited MBH growth with 
meeting the critical condition 



(25) 



where r; csc J> 1 and T v i r is given by equation (1151 . Once 
heated to T cq = rj csc T v i T , the gas will stop feeding the MBH. 
We readily find that the condition (1251 will be met for 



M BH ,crit = 4.6xl0 10 M Q)?CS c f A 

\2Q0kms~V 



,200km s^J L At 
Therefore, for fixed values of »7 CS c, L/Z/Edd and M gas /M we 
expect MflH.crit oc a 4 , similar to the observed Mbh-o" re- 
lation. According to our proposed scenario, once the MBH 
has reached the critical mass, its accretion growth will be 
effectively terminated so that at the present epoch Mbh is 
expected to be only slightly larger than MBH.crit- 

Equally important information is contained in the nor- 
malization of the Mbh-ct relation. In fact, by comparing 
equation 12611 with equation 11611 we find that the observed 
relationship will be established if 

M eas „ .„„ _ n L At 



Ah 



3x 10-Vi- 



M* 



To satisfy the observed Mbh-c relation, the gas-to-stars ra- 
tio is thus required to be relatively low and approximately 
constant for spheroids with different masses at the epoch 
when the MBH reaches its critical mass, although the ob- 
servational uncertainty in the M B h-o" relation leaves some 
room for a weak dependence of M gas /M* on cr. As for the 

3 Quite obviously, this effect will cooperate with the energy input 
from stellar winds and SNII explosions in the forming galaxies, 
however here we show that in principle the radiative feedback 
from the MBH only is already sufficient to do the work. 
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Figure 5. Thick solid lino shows the predicted Mbh~<t relation 
resulting from the requirement that heating of the interstellar 
gas by radiation from the central MBH at the Eddington limit be 
below the level required to drive the gas from the galaxy (T eq < 
Tvi r ). This upper bound on Mbh is based on assuming a constant 
gas fraction of Mgas/M = 0.003 and 77 CS c = 1. The thin solid line 
corresponds to Mgas/M = 0.0015 and r/ e sc = 2. The dashed line 
is the observed Mbh °c <?" 4 relation in the range Mbh = 10 6 — a 
few 1O 9 M0, extrapolated to lower and higher Mbh values from 
equation 1161 . The dotted lines are Mbh <x and A/bh °c °" 5 
laws. 

Eddington ratio, it is reasonable to expect L/L Edd ~0.1-1 
based on hydrodynamic al simulations (CO01 ) and observa- 
tions of quasars (see e.g. iHaiman et al.ll2004l and references 
therein). 

We note that the above calculation was based on a spe- 
cific (r~ 2 ) gas density profile for which £ = L/nr 2 = const. 
This allowed us to avoid specifying the size of the radiative 
heating zone. Being a reasonable assumption for the outer 
regions of the galaxy, the r~ 2 profile is however expected 
to flatten in the central ~0.1f? e region (see, e.g., CO01), so 
that £ will be increasing toward the MBH. Since it can be 
sufficient to overheat the cooling gas in this central zone to 
produce dramatic effects on the subsequent co-evolution of 
the galaxy and its central MBH, our scenario certainly al- 
lows the critical M gas /M* ratio (defined say within i? e ) to 
be several times larger than required by equation 1271 . 

We note that the approximately linear Teq(£) depen- 
dence [see equation £J] was crucial to the above argu- 
ment leading to the MBH.crit oc a 4 result. However, the 
T cq (£) function becomes strongly nonlinear outside the range 
2 x 10 4 K < T cq < 10 7 K, and a more general result can be 
obtained if we consider the exact curve T cq (£) from SOS. In 
Fig. we show the predicted correlation between MBH.crit 
and a for ?7 osc = 1, L/L E dd = 1 and M gaa /M = 3x 10~ 3 , and 
compare it with the observed relationship. We see that the 
Mbh oc <t 4 behavior is expected to break down for Mbh < 



1O 4 M0 and also for Mbh > a few 1O 9 M0. In the same figure 
we demonstrate the effect of the escape parameter ?7 osc . For 
7/ osc = 2 and the gas fraction decreased two- fold, the Mbh-o 
relation is unchanged except that the high-mass cutoff oc- 
curs at a lower mass. Specifically Mbh, cut ~ 3 x 10 9 Mqi^c- 
It is perhaps interesting that the range of masses shown 
in Fig. |S] for which Mbh oc <t 4 is obtained from considera- 
tions of atomic physics (and the observed AGN spectra) 
corresponds closely with the range of masses for which this 
power law provides a good fit to the observations. Explor- 
ing the Mbh-o" relation observationally near 1O 9 M0 would 
be a sensitive test of the importance of radiative feedback. 
We note that other scenarios discussed in the literature also 
predict deviations f rom the power law trend for the most 
massive MBHs fe.g. lMenci et al]l2003h . 



4 A SIMPLE, PHYSICALLY BASED 
TOY-MODEL 

In this Section we address in a more quantitative way the 
MBH-galaxy co-evolution. We adopt a physically-motivated 
one-zone model described in detail in Appendix^ focusing 
in particular on the co-evolution of the galaxy gas budget 
and stellar mass, a key ingredient of the proposed scenario 
(see Sect. 3). Several aspects of this model have been already 
described in Ciotti, Ostriker & Sazonov (2004) and Ostriker 
& Ciotti (2004). 

In particular, source terms for the galaxy gas mass 
(equation IA1I 1 are due to cosmological infall (equation IA2I 1 
and to stellar mass return from the evolving stellar popula- 
tion (equation IA4I 1. while gas is subtracted from the total 
budget by star formation (equation IA3I ). gaseous MBH ac- 
cretion (equation IA6I ). and (possible) galactic winds when 
the thermal energy of the ISM is high enough to escape 
from the galaxy potential well (equation |A7| 1. Energy in- 
put on the galactic gas is due to thermalization of SNII and 
SNIa explosions, to thermalization of red giants winds due to 
the galaxy stellar velocity dispersion (Sects. A3.1 and A3. 2), 
and to the radiative feedback from the accreting MBH (Sect. 
A3. 3). In case of galactic winds, we also consider adiabatic 
cooling of the expanding gas (Sect. A3. 4). Stars are formed 
by gas cooling, while two different mechanisms are consid- 
ered for the MBH growth, namely gaseous accretion and 
coalescence of stellar remnants of massive stars (equation 
IA5I . The galaxy potential well is determined only by the 
dark matter potential well, which is assumed to be unevolv- 
ing with cosmic time. In other words, here we are not con- 
sidering merging between galaxies; finally, we note that in 
the present toy-model the possibility to add a simple recipe 
for chemical evolution is straightforward, so that one could 
also test the model against scaling laws such as the Mg2-<r 
relation. 

We remark here that this kind of approach is not new 
(e.g. iGranato et alJliool iMenci et aljEllll.Ti. However, our 
scenario presents a few but important differences with re- 
spect to the other cases. For example, in their modeliza- 
tion of the AGN f eedback IGranato et alJ <l2004f) (see also 
iMurrav et all 120041) assumed that the main role is played 
by radiation pressure through scattering and absorption by 
dust and scattering in resonance lines, while the energy re- 
lease from the AGN to the ISM is due to thermalization 
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of the kinetic energy of the outflow. In our case, instead, 
the feedback is due to radiative heating. Note that energy 
absorption is typically more efficient in driving winds than 
momentum absorption. 

Obviously, the results of such an approach should not be 
"overinterpreted" : as common in all similar approaches, the 
parameter space of the present model is huge (even though 
several input parameters are nicely constrained by theory 
and/or observations), thus the results of such kind of simu- 
lations should be interpreted more as indications of possible 
evolutionary histories than exact predictions. In particular, 
the toy-model cannot directly test the ability of radiative feed- 
back to produce the right final MBH mass, in fact, this can 
be done only using numerical hydrodynamical simulations. 
This is not surprising, because the toy-model, by construc- 
tion, is a one zone model, and we already know that feedback 
mechanisms are strongly scale-dependent, in the sense that 
central galaxy regions react in a substantially different way 
with respect to the whole system (C097, CO01, Sect. EH . 

As a consequence, one of the key input parameters in 
our scheme is the adopted "quasar duty cycle" /Edd [equa- 
tion (A6)], which essentially introduces a limit, /Edd^Edd, on 
the instantaneous central quasar luminosity that is smaller 
than the Eddington luminosity. In the following, we present 
the results for two different assumptions about /sdd- In the 
first (Sect. 14. H . we distinguish two evolutionary phases: 
a first phase (the "cold phase", in which T coo i/Td yn < 1) 
that would be identified observationally with the Lyman 
Break Galaxies, and a later phase (the "hot phase", in which 
fcooi/Tdyn > 1) that would be identified with normal, local 
ellipticals which contain little gas, have low rates of star 
formation and have a duty cycle (fraction of time during 
which they appear as luminous AGNs) of roughly 0.1%. 
Thus, in equation (A6) the duty cycle factor is known em- 
pirically to be /Edd — 0.01 in the cold phase (where roughly 
a few per cent of Lyman Break Galaxie s show central AGNs, 
ISteidel et ai]l2002l : iLehmer et alJl2004ll . while / Ed d ^ 0.001 
in the hot phase. It is important to note that the last value 
is suggeste d by both numerica l simulations (CO01) and ob- 
servations jHaiman et aljl2004f) . 

A different model is also explored (Sect. l4~2t . in which 
/Edd is kept fixed to 1 over all the simulation. This case will 
adequately demonstrate the effects of radiative feedback on 
the toy-model evolution. 



4.1 Models with reduced /Edd 

With the above remarks in mind, here we present a few 
representative simulations primarily aimed at investigating 
whether there exists a specific evolutionary phase character- 
ized by a gas-to-star mass ratio of the order of the factor in 
equation 12711 . required to obtain the right Mbh-o" relation. 
The time evolution of the quantities shown in Figs.|SJ|S]refers 
to a galaxy reference model characterized by R a = 4kpc 
and a halo (constant) circular velocity of 400kms _1 ; the 
total mass of the gas infall is 10 Mq, and the characteris- 
tic infall time is 2 Gyr. Other relevant simulations param- 
eters are a, = 0.3 in equation 1A3L Pbh,* = 1.5 x 10~ 4 
in equation 1A5I . e = 0.1 in equation 1A17L tjsn = 0.5 in 
equation 1A27L and finally ?7 CSC = 2 in equation (IA7II . The 
initial black hole mass is assumed to be IOMq , and the duty- 
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Figure 6. Reference model. Panel a: mass infall rate (solid 
line) and stellar mass formation rate (dotted line). Panel b: total 
(gaseous plus stellar remnants) MBH accretion rate (solid line), 
Bondi accretion rate (dotted line), Eddington accretion rate (re- 
duced by the duty-cycle factor /Eddi dashed line). Panel c: total 
infall mass (solid line), total stellar galaxy mass (dotted line), to- 
tal galaxy gas mass (short-dashed line); the nearly horizontal line 
is the escaped gas mass. Panel d: total MBH mass (solid line), to- 
tal mass gaseously accreted (dotted line), MBH mass originated 
from stellar remnants (dashed line). 



cycle is fixed according to the prescription of the "cold/hot" 
phases described above. 

While a complete description of the toy-model behav- 
ior for different choices of the input parameters is outside 
the scope of this paper, here we remark that after an ini- 
tial "cold" phase dominated by gas infall, as soon as the gas 
density becomes sufficiently low (Fig. EJ:) , and correspond- 
ingly the cooling time becomes longer than the dynamical 
time (Fig. 01), the gas heating dominates, and the galaxy 
switches to a "hot" solution (Fig. 0. The gas mass/stellar 
mass ratio at that moment (~ 0.003, Fig. 0;) is very near 
with the value inferred in Section[3]from the argument lead- 
ing to the right A/bh-o" relation. Note also how the gas con- 
tent of the present day "galaxy" is in nice agreement with 
observations. It proves that if /Edd is sufficiently small (as 
in the present case) , the accretion is Eddington (rather than 
Bondi) limited (equation [A6]) throughout the MBH-galaxy 
co-evolution, and then our model increases the MBH mass 
in the prescribed way independently of the specific feedback 
mechanism. The obtained final black hole mass is in nice 
agreement with the Magorrian relation (dotted line, Fig. 
7b), but the significance of this result should not be overes- 
timated since in reality the Magorrian relation is set (in our 
scenario) by the MBH feedback. 

An interesting experiment is obtained by reducing the 
circular halo velocity and the infall mass in the reference 
model: in these cases galactic winds are favoured. In other 
words, small galaxies lose their gas content easily, in accor- 
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Figure 7. Reference model. Panel a: time evolution of the Bondi 
radius. Panel b: logarithm (base 10) of the ratio between gas mass 
to stellar mass (solid line) and of the ration between MBH mass 
to stellar mass ("Magorrian relation", dotted line). Panel c: gas 
density at the Bondi radius (solid line) and mean gas density 
(dotted line). Panel d: logarithm (base 10) of the cooling time 
(solid line) and heating time (dotted line) measured in terms of 
the dynamical tyme. 




t/Gyr 

Figure 8. Time evolution of the model gas temperature (solid 
line). The model virial temperature is represented by the dotted 
line, while the dashed line represent the "escape" temperature 
(here assumed 2T v j r ). 



dance with the Mg2-er and the Faber-Jackson relations and 
with the hydrodynamical simulations of COOl. Remarkably, 
the transition to the hot phase of these models happens for 
Mgas/M* ~ 0.01, similarly to the behavior seen in the case 
of the more massive spheroid in the reference model. In other 
models we have verified the role of the AGN photoionization 
by artificially excluding the sources of stellar heating in the 
reference model: in this case the galaxy gas temperature (af- 
ter the initial cold phase) remains sub-virial (i.e. below the 
horizontal dotted line in Fig. 8), without however cooling 
down to the imposed lower temperature limit of 10 4 K, as 
instead happens if we also exclude the MBH feedback. This 
sort of cooperation between AGN feedback and stellar en- 
ergy injection, i.e., the fact that substantial galactic winds in 
general are due to stellar heating, and are reinforced by the 
presence of the central A GN, was already found in numerical 
simulations (COOl). 

An important and apparently robust conclusion that 
can be drawn from these simulations is that stellar heating 
inevitably leads to a transition from cold to hot solution 
when the gas-to-star mass ratio drops to of order 1 per cent 
or somewhat less. Now, since a gas fraction of this order is 
required for the radiative feedback from the central MBH to 
limit its growth at the mass obeying the observed Mbh-o" 
relation (see Sect. El, it is tempting to suggest that the MBH 
reaches its critical mass, determined by radiative feedback, 
approximately at the epoch of transition from cold to hot 
galaxy phase. 

We emphasize again that in the above models the radia- 
tive output from the MBH does not strongly influence the 
interstellar gas, since /Edd is assumed to be low and thus the 
central luminosity is forced to be less than /Edd^Edd -C ^Edd 
at any instant. In other words, these models allow persistent 
moderate AGN activity but forbid quasar-type outbursts. 



4.2 A model illustrating radiative feedback 

In order to illustrate the effects of MBH feedback in the 
context of the toy model, we now present a model in all 
similar to reference model (Sect. POll . but in which /Edd = 
1 over all the evolution. In this model we also reduce the 
characteristic infall time from 2 Gyr to 0.2 Gyr, and we set 
Pbh,* = in equation 1A5L i.e., we are neglecting black 
hole growth due to accretion of stellar remnants. 

Several differences are apparent with respect to the evo- 
lution of the reference model and its variants. In fact, the 
final black hole mass is now larger, the final amount of gas 
in the galaxy is significantly lower, and the galaxy is found 
in a permanent wind state (excluding the very short initial 
cold phase). Also, from Fig. 9b is apparent how the MBH 
grows by Bondi accretion instead of being Eddington lim- 
ited (again excluding the initial phase). The main difference, 
however, is that now a very strong gas outflow is caused by 
the MBH feedback rather than by energy input from the 
evolving stellar population. 

Most importantly, this model supports the argument 
presented in Section |3] in that quasar radiative heating ter- 
minates the MBH growth at a mass proportional to the gas 
fraction at the critical epoch, i.e. when the MBH feedback 
becomes important. The fact that in this particular simu- 
lation the final black hole mass turns out to be about 2 
orders of magnitude higher than the "Magorrian" mass is 
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Figure 9. Model with /Edd = 1 an d 7"in,f = 0.2 Gyr. See Caption Figure 10. Same model as in Fig. 9. See Caption of Fig. 7 for a 

of Fig. 6 for a description of the various curves. Note also how in description of teh various curves, 

the hot phase MBH accretion is Bondi dominated. 



due to the fact that the radiation from the central quasar 
overheats the gas when its mass fraction is ~ 10 per cent 
(see Fig. 10b), much higher than a fraction of 1 per cent. 
This happens in this model because the Eddington limited 
growth of the black hole allowed it to reach a very high mass 
while the galaxy gas content has not yet been reduced to low 
levels. 

An obvious problem with the models described in Sec- 
tions 14.11 and 14.21 is that the parametrization of the MBH 
accretion rate and luminosity in terms of the time aver- 
aged /Edd factor is a very poor substitute for considering the 
real, time dependent problem. Nevertheless, we believe that 
the models assuming a small /Edd and the model in which 
/Edd = 1 nicely complement each other and provide impor- 
tant clues to the real picture of MBH-galaxy co-evolution. 
Specifically, the first kind of models demonstrate that there 
is indeed a well defined phase in galactic evolution when 
the gas-to-stars mass ratio is of the order of 1 per cent - 
it corresponds to the transition from cold to hot solution. 
On the other hand, the /Edd = 1 model demonstrates how 
the MBH radiative feedback can be efficient if the central 
quasar switches on at a near Eddington luminosity. Taken 
together, these models suggest that a major quasar outburst 
occuring approximately at the beginning of the hot galaxy 
phase, when M gas /M* ~ 0.01 or somewhat less, can lead to 
a significant degassing of the galaxy and termination of the 
MBH growth at the mass required by the Mbh-o" relation. 



5 DISCUSSION AND CONCLUSIONS 

In this paper we explore the possible role of AGN radia- 
tive feedback in the co-evolution of MBHs at the center of 
spheroids and their stellar and gaseous components. The 
feedback is due to a combination of photoionization and 
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Figure 11. Temperature evolution of the model in Figs. 9 and 
10. See Caption of Fig. 8 for a description of the various curves. 
Note how the galaxy is found in a permanent wind phase after 
~ 1 Gyr. 



Compton heating. In our calculations we adopt (from the 
work of SOS) a broad-band spectral energy distribution 
that corresponds closely with the radiative output of typ- 
ical quasars. 

We first investigate (Sect. ETTI . based on energetic con- 
siderations, whether the radiative output of the central 
MBH can significantly affect a Bondi-type accretion flow in 
the central 1-100 pc of a galaxy. In agreement with previous 
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studies, we find that if the accretion is radiatively efficient 
and proceeds at a subcriticai rate, as expected of quasars, 
the mass inflow rate and the MBH iuminosity are expected 
to osciilate on time scales ~ 10 4 -10 5 yr as a result of gas pre- 
heating, unless the temperature of the supplied interstellar 
gas approaches the Compton temperature of the radiation 
field (2 x 10 r K). We also note that if a central accretion 
disk extends beyond ~ 0.001-0.01 pc from the MBH, then 
preheating of the external, hot flow may still lead to varia- 
tions in the MBH accretion rate, but on longer time scales 
(up to millions of years), determined by the characteristic 
drift time in the disk. 

Using similar energetic considerations, we then assess 
the effect of the MBH radiative feedback on the main vol- 
ume of the host galaxy (Sect. 12.21 . We demonstrate that 
the radiative output from a MBH growing by radiatively 
efficient accretion can unbind the ambient interstellar gas, 
provided it is tenuous enough, out to a maximum distance 
of a few kpc, which would typically amount to a significant 
fraction of the effective radius for giant ellipticals and to a 
few effective radii for smaller spheroids. 

We then discuss (Sect. [3J a possible origin of the ob- 
served Mbh-c relation based on the hypothesis that the ra- 
diative output from the growing MBH will eventually heat 
the ambient interstellar gas above the virial temperature 
and expel most of it, limiting both future growth of the 
MBH and future star formation to low levels. We demon- 
strate that if the gas-to-stars mass ratio drops to <1 per 
cent in the central regions of spheroids at a certain stage of 
their evolution, then the radiative feedback from the cen- 
tral MBH switched on as a bright quasar will terminate its 
growth at the mass obeying the observed Mbh-o" relation. 
Furthermore, we predict that the observed power law should 
break down for black hole masses below 10 4 M© and above 
a few 10 9 Mq. These considerations however leave open the 
question as to why the gas fraction should be of the above 
order when the MBH reaches the critical mass. 

In order to obtain a better insight into the co-evolution 
of gas, stars and the central MBH over the Hubble time we 
explore a simple but physically based one-zone toy model 
(Sect.QJ. A robust result obtained from the simulations is 
that when the gas mass fraction has been reduced by star 
formation to of order 1 per cent or somewhat less, a transi- 
tion from cold to hot solution takes place, both in the case 
of massive protogalaxies that would evolve into the present 
day giant ellipticals and in the case of smaller spheroids. 
This transition occurs primarily as a result of gas heating 
due to the evolving galaxy stellar component in its various 
forms (stellar winds thermalization, SNII and SNIa explo- 
sions) , while the AGN heating is also expected to contribute 
at some level. 

The near coincidence of the gas fraction correspond- 
ing to the beginning of the hot galactic phase with that 
[equation 1271 1 required by our argument leading to the cor- 
rect Mbh-o" relation offers the possibility of the following 
evolutionary scenario. At the early stages of galaxy evolu- 
tion when the protogalactic gas is dense and cold, active 
star formation is accompanied by the growth of a central 
black hole. The black hole however is not massive enough to 
produce a strong heating effect on the ambient, dense gas, 
even during episodes when it shines near at the Eddington 
limit. This cold phase would be identified observationally 



with the Lyman Break Galaxies and bright submillimeter 
galaxies, which are characterized by high star formation rate 
and moderate AGN activity. The cold phase ends when the 
gas-to-stars mass ratio has been reduced to ~ 0.01, when 
the energy input from the evolving stellar population and 
possibly from the central MBH heat the gas to a sub-virial 
temperature. The MBH continues to grow actively during 
this transitional epoch (perhaps more actively than before) 
because there are still sufficient supplies of gas for accre- 
tion, and soon reaches the critical mass (obeying the Mbh~ 
a relation), when the MBH radiative output causes a major 
gas outflow. This phase would be identified with the ma- 
jor quasar epoch. The subsequent evolution is passive and 
characterized by very low AGN activity and a duty cycle re- 
duced by a factor of ten to 0.001, and this late phase would 
be identified with the present day elliptical galaxies. 

This scenario is admittedly tentative and not the only 
possible one, and should be verified by observations and 
more detailed computations. In particular, it should be eluci- 
dated as to why the MBH does not grow to an excessive mass 
during the cold galactic phase, a key assumption implicitly 
made above. Possibl y this does not happen , as also suggested 
in other works fe.g. lArchibald et alJl20o'3) . just because the 
stellar spheroid is formed too rapidly even compared to an 
Eddington limited growth of the MBH from an initial stellar 
mass up to 1O 8 -1O 9 M0, which takes <JlGyr. Alternatively, 
early MBH growth may be characterized by some duty cycle 
determined by the poorly understood physics of gas supply 
and accretion on to protogalaxies. 

A proper investigation of the importance of radiative 
heating on the MBH-galaxy co-evolution, based on high spa- 
tial hydrodynamical numerical simulations and adopting the 
specification of the input spectrum and atomic physics from 
this work, is now in progress (Ciotti & Ostriker, in prepara- 
tion). 
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APPENDIX A: A TOY MODEL FOR THE MBH 
GROWTH 

Al The equations 

The equations describing the evolution of the physical quan- 
tities considered in the present toy-model are 



M 



M, 



inf 



M* + M r 



M B 1 



M e . 



(Al) 



M lnf = — — exp I 



M„ 



ce,M aas 



ma,x(r dy 



M re 



(A2) 
(A3) 
(A4) 



M rec = / M+{t')W*{t - t') dt', 
Jo 

where Mt is the first term on the r.h.s. of equation l|A30 

M BH = Mbh.ucc + /W»M+, (A5) 
where 

M B H,acc = min(/ E ddM E dd, Mb), (A6) 
and finally 



Me 



(A7) 



0, T < ri esc T vil 



A2 Input physics 

A2.1 Gas equilibrium distribution 

The code is started by assigning the dark-matter halo circu- 
lar velocity v c under the assumption of a (singular) isother- 
mal distribution [see equation 123H 1. and the quantity Td y n 
entering equation <A3t is defined as 



T~dyn 



2ttR c 



(A8) 



where R e is a characteristic scale-length that could be iden- 
tified with the effective radius or the half-mass radius of the 
gas and galaxy stellar distribution. For example, the radial 
trend of the gas density is obtained by arbitrarily impos- 
ing that all the gas mass is contained within 2R a , and it is 
distributed as the dark matter halo: 



P = 



pe_ ( Rq\ 
3 V r ) 



(A9) 



The assumption above should not taken too literally, in the 
sense that it is only a simple way to obtain a representative 
value of the mean gas density p e : different choices would 
lead to different values of p e for an assigned total gas mass. 
In any case, as a consequence of our assumption, the mean 
gas density within R B is given by 



3 fl^Igas 

8nR% 



(A10) 



The (mass- weighted) equilibrium gas temperature T v i r could 
be obtained from the hydrostatic equilibrium and Jeans 
equations. However it is well known that a gas distribution 
as that in equation &A91 but untruncated (and so character- 
ized by an infinite mass) would have the (constant) temper- 
ature 



fj,m p Vc 



/, 2k • ' Alli 

where we recall that for the singular isothermal sphere a 2 — 
v 2 /2, where a is the (1-dimensional) halo velocity dispersion, 
and pi — 0.61. For simplicity we use equation 1A1H as an 
estimate of the equilibrium gas temperature associated with 
the density distribution in equation (IA9ft . 
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A2.2 Gas cooling time 

The gas mean cooling time within R a is defined as r e0 oi 
E/Ec, where 



E = 

2/im p 

is the gas internal energy per unit volume, 



(A12) 



A(T), (A13) 



^ = n e n p A(T) = !jA(T)=^ 
where 

MT) = 2A8 T oslV S + 2 -706 x lO^T 2976 ergcm 3 /s(A14) 



fsee lMathews fc Bregmanlll978l CO01). We thus obtain 
6fim p k T 



Tcool 



p e A(T) • 



(A15) 



A2.3 Stellar physics 

The gas recycled by the evolving stellar population is given 
by equation 1A4I . The modulating (normalized) kernel is 
given by 



W,(t) = R*x 



To 



t + TQ 



(A16) 



where we adopt R t = 0.3, 5» = 1.36, to = 10 8 yr; thus, a 
fraction R, of the stellar mass produced at any time step is 
recycled at the end to the galaxy interstellar medium. The 
function above is a fairly acceptable fit of the time dependent 
mass return rate of a simple, pas s ively evolving stellar pop- 
ulation, adopted in ICiotti et alJ l|l99lfr (hereafter CDPR), 
C097, 01. 



A 2. 4 MBH accretion physics 

Concerning equations 1A5I and 1A6L we define 

-^Edd 



M] 



(A17) 



where 0.001 < e < 0.1 and L-^dd is given in equation 
The Bondi accretion rate is given by 



M B = 4nR B p B c B , 

where _Rb is given by equation 10 and 

'dp\ /,-:/' 



d PJisot I™*' 



(A18) 



(A19) 



An estimate of the gas density at Rb, Pb, is obtained from 
equation IA9t evaluated at r = Rb 4 - 



4 Note that the Bondi accretion rate given by equations above 
becomes irrelevant when the gas is colder than the virial tem- 
perature, because the Bondi radius is then larger than the MBH 
radius of influence. More appropriate in this case would be to 
calculate the accretion rate for the conditions at the influence ra- 
dius. However, for the adopted gas density profile the accretion 
rate will in fact remain unchanged. 



A3 Feedback and gas escape 

In the following two sections we describe two different kinds 
of heating that act on the galactic gas, plus the gas cooling: 
note that we are not counting cooling twice, because the 
cooling function in equation <A13t is used to determine the 
gas cooling time only. Thus, we determine the gas temper- 
ature at each time-step by integrating the equation for the 
internal energy per unit volume 



E = E 



H,SN 



+ Eh,w + Eh, 



-E c + E ad + 



t M inf \v 



2 

esc 



" 16^1 ' (A20) 

where A is a dimensionless parameter ranging from 0.25 (cor- 
responding to virial velocity) to 1; E a d describes the adia- 
batic cooling in case of gas escaping (Sect. A3. 4), while Eh,w 
describes the energy input due to the thermalization of red- 
giants wind (Sect. A3. 2). As a reference value we use the 
quantities within the half mass radius (here identified with 
R e ) and, accordingly, the gas internal energy is given by 
equation (IA 121 . while only half of the gas accretion/outflow 
kinetic energy is considered. In particular, in the integration 
of equation above we adopt a global fitting formula for the 
quantity Eh.agn — E c (Sect. A3.3). From equations <A"T2t 
and dA^Ot 

2pm p E(t + At) 



T(t + At) 



(A21) 



3kp e (t + At) 

At each time-step a check on the attained gas temperature 
is performed, and we force the gas temperature to remain 
above the minimum value T m i n = 10 4 K. Energy (temper- 
ature) integration is actually the last time-step in the inte- 
gration cycle. When a new interation starts, we first check if 
the new determined temperature is higher than the escape 
temperature 77escTvj r . If the condition is not satisfied (i.e., 
T < ryoscTvir), gas escape is suppressed. On the contrary, if 
T > fjescTvir, we adopt 

Tesc = ~ (A22) 

Cs 

in equation 1A7I . 

Note that an elementary application of the Virial The- 
orem, coupled with the definition of escape velocity and un- 
der the assumption that the gas is distributed as the dark 
matter halo, would predict that its mass weighted escape 
temperature correspond to rfcsc = 4. However, this is cer- 
tainly an upper limit, and we expect that feedback effects 
will be important at much lower temperatures, say ~ 2T v i r 
or even less. 



A3. 1 SNII, SNIa, OB and red giants wind heating 

We describe here the first two contributions to the gas heat- 
ing in equation 1A201 . The supernova heating Eh,sn = 
En + Eia is given by the sum of SNII and SNIa energy 
injection. Due to the very short life-time of massive stars, 
the SNII heating is assumed istantaneous. 

We now compute the expected number of SNII events 
from a new stellar mass added to the galaxy according to the 
first term in the r.h.s of equation <A3t . Let M* be the total 
stellar mass assembled in the unit time, without considering 
the successive mass losses. We adopt a Salpeter IMF distri- 
bution, with the standard low mass cut-off M*,i = O.IMq, 
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while for simplicity we assume an infinite value for the upper 
mass: 

* = = M+(x - VjM*- 1 M- (1+x) , (A23) 

with x = 1.35. The number of SNII events per unit time is 
just given by the number of stars with M > Mri = 8Mq : 



Nr 



M T1 



V xj M,a \M n J 



(A24) 



The total energy released by the SNII per unit time is then 
given by NuEsn (where Esn = 10 51 erg/s), while the ac- 
tual gas heating is tjsnNiiEsn, where rjsiy is an efficiency 
factor. Note that the value of the efficiency t]sn is highly 
uncertain, depending on the particular status of the ISM. In 
addition, in our simple formalism in the factor t]sn > use 
also consider the energy injected by stellar winds by the mas- 
sive OB stars. In the computed models we adopt rjs jv = 0.5. 

The treatment of SNIa heating is more complicate, be- 
ing caracterized by a delayed explosion with respect to the 
star formation episode. We adopt here the same description 
as in CDPR, C097, COOl. Thus, if N Ia dt' is the total num- 
ber of SNIa that will explode at t > t' as a consequence of 
the star formation episode at t', and Wi a {t — t') the normal- 
ized explosion rate, one obtains 



Eia = r/sNEst 



N Ia {t')Wi a {t-t')dt', 



(A25) 



where for simplicity we adopted the same efficiency as for 
SNII. Following CDPR 

Sla - 1 



Wla{t) = 



Tla 



(A26) 



Tla \t + Tla 

where Sj a — 1.4 — 1.6 and we arbitrarily assume n a — 0.5 
Gyrs. 

We do not enter here in the complicate problem of the 
determination of the present rate and the total number of 
SNIa explosions in a given stellar population: in the present 
treatment we just assume that for each of the star forma- 
tion episodes Ni a = $i a Nn. The actual value of i?/ a is fixed 
by (arbitrarily) requiring that the standard solar propor- 
tion is maintained, i.e., that SNIa provide 3/4 of the iron 
in the ISM, while t he remaining 1/4 is provided by SNII 
jRenzini et al.lll993ft . Because each SNIa produces ~ 10 
times more Iron that a SNII, then i?/ a ~ 0.3. 

Thus, the gas heating per unit volume within R e due to 
SNII, stellar winds from OB stars, and SNIa, is given by 



E. 



3t]snE S n(Nii + Nia 



(A27) 



8tt.RI 

Finally, thermalization of stellar winds emitted by red-giants 
to the stellar velocity dispersion "temperature" is an impor- 
tant contribution to the global energy budget of the ISM in 
early type galaxies (see, e.g., CDPR). Here we describe this 
term as 

Eh - w -^w (A28) 

i.e., we evaluate the energy input per unit volume inside _R C . 

Thus, the characteristic heating time due to the global 
energy input of the evolving stellar population is 

3kp e T/2fim p 



Theat,SN,w 



E, 



(A29) 



A3. 2 Adiabatic cooling 

If the gas escape condition is satisfied, the escaping gas not 
only carries out kinetic energy, but also decreases its in- 
ternal energy in order to expand. We evaluate the cooling 
due to gas expansion by using the First Law of Thermodi- 
namics. In particular, the mechanical work of expansion is 
given by — pdV. During the expansion the gas mass is con- 
served, so that V = Mg as /p and dV = —M gaa dp/p 2 . The 
variation of internal energy (per unit volume, as required in 
energy equation above) due to gas expansion is then given 
by pM gas dp/(Vp 2 ) = pdlnp. Then, from p oc r(t) -3 , and by 
assuming that the expansion velocity is given by the sound 
velocity, we obtain 

Ead = -e£ = -2—, (A30) 

^Cc Te.sc 

where we used the identity E = 3p/2 and as order of magni- 
tude r ~ 2R a . The equation above is integrated with time- 
splitting. No adiabatic heating is considered in case of infall, 
being considered as accretion of clouds. 



A3. 3 AGN radiative feedback 

Using XSTAR <Kallmanl i2002'l we compute the net volume 
heating rate E of a cosmic plasma exposed to radiation char- 
acterized by the average quasar SED adopted from SOS. 
Photoionization equilibrium is assumed. E depends on gas 
temperature T and on the ionization parameter f , which can 
be determined from equations and &A91 as follows: 

^^^ = JM' (A31) 

where L — €C 2 Mbh is calculated at each time-step from 
equation ijABfl . Note that from equation l|A9|l £ is indepen- 
dent of radius (as far as the gas is optically thin) . 

The following furmula provides a good approximation 
to E (all quantities are expressed in cgs system): 

E = n 2 (S 1 + S 2 + S 3 ), 



where n is the H nucleus density (in number) 

5 1 = -3.8 x io _27 Vr 

are the bremsstrahlung losses, 

5 2 = 4.1 x 10~ 35 (1.9 x 10 7 — T) £ 
is Compton heating (cooling), and 

Ss = 10 -23« + &(C/eo) C 



(A32) 
(A33) 
(A34) 

(A35) 



i + (£/£o) c 

is the sum of photoionization heating, line and recombina- 
tion continuum cooling. 
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e 25(logT-4.35) 2 e 5.5(logT-5.2) 2 
1.7 X 10 4 



3.6(logT-6.5) 2 ' 



yo.' 



C= 1.1 



1.1 



oT/1.8 10 5 



+ 



4 x 10 1 



T 4 



(A36) 
(A37) 
(A38) 



and finally 



1.5/VT+ 1.5 x W 12 /VW 
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4 x 10 

y2 



10 r 



80 



1 + 



e (T-10 4 )/1.5 10 3 



(A39) 



The formula above is applicable in the temperature range 
10 4 < T < 3 x 10 7 independently of the £ value, except for 
T < 2 x 10 4 when it breaks down at £ < 0.01 (hydrogen 
becomes neutral), but then Si, S2, S3 — > 0. 

For T > 3 x 10 7 K one can use the approximation E = 



n 2 (Si + S* 2 ). 



